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We study the higher-order correlation functions of covariant families of observ- 
ables associated with random Schrodinger operators on the lattice in the strong 
disorder regime. We prove that if the distribution of the random variables has a 
density analytic in a strip about the real axis, then these correlation functions are 
analytic functions of the energy outside of the planes corresponding to coincident 
energies. In particular, this implies the analyticity of the density of states, and of 
the current-current correlation function outside of the diagonal. Consequently, this 
proves that the current-current correlation function has an analytic density outside 
of the diagonal at strong disorder. 
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1 Correlation Functions 

The study of the higher-order correlation functions for random Schrodinger operators is 
essential for an understanding of the transport properties of the system. In this paper, 
we study the correlation functions for covariant families of observables associated with 
random Anderson-type Schrodinger operators on d- dimensional lattices Z d in the strong 
disorder regime. The Anderson model is given by the following random Hamiltonian 
acting on £ 2 (Z d ) 

H u ip(x) — A Yl ^(V) + V u (x)ip(x), for ^ef(Z d ). (1) 

y;\y-x\=l 

Here A is a small real parameter providing a measure of the disorder of the system. The 
random potential V u is determined by a family of independent, identically distributed 
(iid) random variables to = (Vu(x)) xeZd with a common distribution given by a density 
g(V)dV. We assume that the density g is analytic in a strip about the real axis. We 
write P for the probability measure generated by g as an infinite product measure on the 
probability space fl Precise assumptions and formulation of the results are given in the 
next section. 

In order to describe our results, let us suppose that we have a family of covariant 
observables A Wi i, for i = 1, . . . , N. A basic example is the j t/l -component of the velocity 
operator Vj = i[H w ,Xj], that is independent of uj. The resolvent Rw{z) = (H u — z)' 1 
for H w has matrix elements G w (x, y;z) = (x\R UJ (z)\y) giving the Green function. The 
associated spectral density operator p u (E) is given by lim e ^ $S{H U — E — ie)^ 1 . The Ap- 
point correlation function Kpj, associated with a covariant family {A^j \j = 1, . . . , N}, 
is given by 

K N (E U E 2 ,...,E N )= [ dF(u)(0\p„(E 1 )A,, tl p„(E 2 )A w>2 ---p„(E N )A w>N \0). (2) 

Jn 

We are interested in the behavior of Kn as a function of the energies Ej, j = 1, . . . , N. 
It is not a priori clear that Kn in (J2J is well-defined. We will prove, in fact, that 
Kjy(Ei, . . . , En) is real analytic in a region of ~R N away from the planes where Ej = Ei, 
for i 7^ j. 
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As defined in (j2J), the first correlation function Ki(E), traditionally noted as p{E), is 
the density of states (DOS) of the family H^. This function has been extensively studied, 
especially in one-dimension, and we refer the reader to the monographs [SJ El for results 
on the regularity of the DOS for one-dimensional lattice models. In this article, we are 
especially concerned with lattice models in d > 2 for which there are very few analyticity 
results. Constantinescu, Frohlich, and Spencer [TUj studied the DOS for lattice models 
using the random walk expansion described in section 3. First, these authors proved 
two results independent of the (nonzero) disorder. They proved that if g is analytic 
in a strip of width 5 > 4d, then p(E) is analytic in a small strip around \dtE\ >> 1. 
They also proved that if g is Gaussian, then p(E) is analytic for \dtE\ » 1, in a region 
\$sE\ < (l/y2)\3t,E\. Furthermore, the DOS decays like a Gaussian as \E\ — > oo along the 
real axis. Note that both of these results are large energy results. Secondly, they proved 
that for g Gaussian and large disorder, the DOS p(E) is analytic in a neighborhood of the 
real axis R. Later, Bovier, Campanino, Klein, and Perez [7] obtained stronger smoothness 
results on the DOS using a supersymmetric representation of the Green's function. These 
authors (actually, attributed to Klein and Perez, unpublished, in [7. ) prove two classes 
of results of interest here. First, they prove that if the characteristic function h of the 
probability measure satisfies (1 + t) {d+k) h(t) E L\R), then N(E) G C (fc+1) (R), for all 
disorder sufficiently large (including k = oo). Secondly, if e at h(t) is bounded for some 
a > 0, then for any < a± < a, there is a constant Ai > so that for any < A < Ai, 
the IDS N(E) is analytic in a strip \QE\ < aci- Our results, restricted to the case N = 1, 
are essentially the same as this last mentioned result on the analyticity of the DOS in a 
strip about the real axis provided the density has a continuation. Our method of proof is 
completely different and generalizes to any higher-order correlation functions. 

The higher-order correlation functions have not been studied systematically, although 
L. Pastur [T^], and one of the authors [HIE], have long emphasized their importance in 
transport problems. The second-order correlation functions K2(E±, E2), for the choices 
— Vi, are called the current- current correlation functions. These functions (actually 
measures) play an important role in the theory of conductivity. Until now, it was not 
known if the measure has a density. Our result on the current-current correlation function 
is the first proving the existence of a density for this measure away from the diagonal. 
We will comment further on this in section 2. The second-order correlation function 
has been studied by Pastur and Figotin for a one- dimensional quasi-periodic Schrodinger 
operator [TT|. The strong disorder expansion of the second moment suggested by Pastur 
was developed in Kirsch, Lenoble, Pastur providing insight into the behavior of the 
function. In general, very little is known about this function and one of our results is that 
this function is analytic away from the diagonal E\ — E- 2 . The DC conductivity requires 
control of the two-point correlation function on the diagonal Ei — E 2 . Although we do 
not achieve this result in this paper, we give the first results in the study of the correlation 
functions. 

The higher-order correlation functions Kn correspond to not-necessarily-positive, bounded, 
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Radon measures K N (dEi, . . . , dE N ). Indeed, for any N, if /i, • ' ' > /jv are smooth func- 
tions with compact support on K then 

%{V jl Hf 1 (H)V j2 Hf 2 (H)---V jN Hf N (H)} = f K N (dE 1 , . . . ,dE N ) frfa) ■ ■ ■ f N (E N ) , 

(3) 

where jk E {l,...,d}. The left-hand side of (JSj) defines tempered distribution in N 
variables. Moreover, the left-hand side is multilinear and bounded with respect to the 
sup norm of the /j's since the operator VjH is bounded for lattice models. Therefore 
belongs to the dual space of (^(R)®^ — Co(R N ), namely it is a bounded Radon measure 
on Mr. Except for the cases N = 1 and N = 2, corresponding to the density of states 
and the current- current correlation function, respectively, this measure is not necessarily 
positive. 

Control of higher-order correlation functions seems to be necessary, for example, in 
order to control growth in a dissipative model describing the Mott variable-range hopping. 
The higher-order moments of the position operator Rj are also controlled by the higher- 
order correlation functions. To see this, we note that 

Rj(t) - Rj{0) = -% [ ds U,[H, Rj]U* = -i [ ds U s VjH U*, (4) 
Jo Jo 

where U s = e~ lsH and V jH = —i[H,Rj] (see ©)• Consequently, there is a function 
F t (Ei, . . . , E 2n ), depending on t, so that the (2n) t/l -moment of the position operator 
Rj{t) — Rj(0), localized to an energy interval A by projector may be written as 

MPA\R,(t)-RM\ 2np A) = I K 2n (dE 1 ,...,dE 2n )F t (E 1 ,...,E 2n ), (5) 

il 2 ("-')xA 

where 

r.(g 1 ....,^)=f r'TV g, - i)1 - 

Our results on the current-current correlation functions, and the higher-order correlation 
functions, are the first results proving the smoothness of these correlation functions outside 
of the coincident planes. 

1.1 Contents of the Paper 

We state the main hypotheses on the models and the main results in section 2. We 
introduce the random walk expansion of the resolvent in section 3. In section 4, we give 
some basic estimates on Cauchy-type integrals. To illustrate our method, we apply these 
to prove the analyticity of the DOS in section 5. We extend these estimates to multiple 
Cauchy-type integrals in section 6. Analyticity of the boundary values of the Cauchy- 
type integrals is proved in section 7. The main result on the analyticity of the correlation 
functions is proved in sections 8. 
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2 The Models, Hypotheses, and the Main Result 



We now provide precise hypotheses and formulation of our results. Throughout this paper, 
we will assume that the single-site probability density g satisfies the following properties: 

1. Analyticity: g can be continued as a holomorphic function in the strip B r = {z G 

C ; \Qz\ < r} 

2. Boundedness: ||g|| r := sup| w | <r J R dv \g{y + iw)\ < oo. 

The space of functions fulfilling these conditions will be denoted by 7i r and it will be 
endowed with the norm || • || r . Note that if < r' < r, then the norm satisfies ||g|| r ' < \\g\\ r , 
for g G TC r . In addition, it is required that g define a probability distribution with zero 
mean and a finite second moment, namely g(v) > 0, for nel, and 



dv g(v) = 1, / dv g(v)v = 0, / dv g(v)v = M g < oo. 

Jr Jr 

The measure space in which u lives will be denoted by Q and the probability measure 

defined by the infinite product of g{V)dV will be denoted P. Obviously Z d acts on fl by 

bimeasurable bijective maps and P is Z d -invariant and ergodic. A covariant operator is a 

family A = (A^g^ of bounded operators on S) = i 2 (Z d ) such that 

1. The map to G VL i— > A w G S(^) is strongly P-measurable 1 . 

2. If T(a) denotes the unitary operator of translation by a then T(a)A UJ T(a)^ 1 = A^, 
if T a is the action of a G 7L d on Q. 

The Anderson Hamiltonian H^, defined in is an example of a family of covariant 
operators whenever the distribution of the random variables V^(x) has compact support. 
However, due to the analyticity condition on g, the Hamiltonians are almost surely 
unbounded and selfadjoint. It is required then that the field of resolvents u G Q h- > 
(zl — Hu)^ 1 G 23(f)) be measurable. The set of covariant operators is a von Neumann 
algebra 2 with the pointwise algebraic operations (sum, product, adjoint) and with the 
norm \\A\\ = P-essup wen || Aj|| 0- An unbounded, selfadjoint, covariant operator, such as 
H^, is said to be affiliated to this algebra if all its spectral projections for bounded Borel 
subsets of K belong to this algebra. There is a canonical trace on this algebra given by 

%{A) = [ d¥(uj)(0\AjO) = lim -|-Tr A (AJ , P a. e. 
Jn M% d |A| 



1 Note that the Borcl structure of 23 (io) is the same for the weak and strong topologies since ij is 
separable. 

2 A von Neumann algebra is a C*-algebra with a predual |18j . 
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This von Neumann algebra is denoted by L°°(7ip) 3 . Let R = (Ri, • • - , Rd) be the position 
operator acting on as a multiplication operator. Then a derivation V on the algebra is 
defined by 

V = (d ir -- ,d d ), with (d^ = i[R^AJ[, (6) 

on the dense subalgebra of covariant operators so that the operators (d^A)^ are bounded 
almost surely. The Integrated Density of States (IDS) has been proved to satisfy the 
Shubin-Pastur formula (see "2(11 IT7| IT| 12*]) 

J\f{E) = lim —^{eigenvalues of E} = %>{x{H < E)) a.e. u, 

AtZ d |A| 

where \{H < E) denotes the spectral projection of the covariant Hamiltonian on the 
interval (— oo, E]. The density of states (DOS) is the Lebesgue-Stieljes measure dhf . 

Transport properties are expressed through the current- current correlation function 
defined as the positive matrix valued measure on M 2 given by 

%{f{H)d u Hg{H)d v ,H) = [ m uy (dE,dE') fXE)g(E'), (7) 

JRxM. 

for f,g G C (1R). Since the electric current operator is given by J = (e/JfjVH, the 
measure m is formally given by m u y(dE,dE') = (E\J u \E'){E'\J u >\E)dEdE' (up to a 
multiplicative constant) in terms of the matrix elements of the current in the eigenbasis 
of the Hamiltonian. 

One of the main problems for transport in such system is to show that this measure 
is absolutely continuous with smooth density. We prove a first result in this direction. 
Namely, we prove that the density exists and is analytic outside of the diagonal at all 
energies in the strong disorder regime. More precisely, let A = {Ai, . . . , A^} be a finite 
family of covariant operators as defined above, such as the velocity operators. The iV-point 
Green function associated with A is defined by 

G A (z) = % ({H - z l )- 1 A 1 ■■■{H- z N )- l A N ) , z = {z u ■ ■ ■ , z N ) G (C \ R) N . (8) 

Note that the spectral function is a linear combination of two Green functions with 
complex-conjugate imaginary parts. The DOS is recovered by taking N — 1, A\ — 1, 
and taking the imaginary part. The two-point correlation function is recovered by taking 
iV = 2, A\ = A2 = d u H, expressing the spectral family in terms of the Green's functions 
via Stone's formula, and by considering the discontinuity of Ga(z) along Z{ = Ei G KL 
The discontinuities are expressed in terms of the boundary values defined as follows: let 

3 For an intrinsic definition without the help of f), see 0E] 
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o_ = (<7i, • • • , <tat) G {+, — } N , and let e = (ei, • • • , ejv) G . Then, with the notation 
a ■ e — (cr^i, • • • , ctn£n), we define the boundary values of (JS]) as 

G%(E) = \imG A (E + ta-e). 

The iV-point correlation functions Kn defined in (J2J) can be expressed as a linear combi- 
nation of these boundary values. Our main result is the following theorem. 

Theorem 1 Let be the Anderson Hamiltonian defined in (QJ) with a distribution g G 
7i r . Then, given N G N, there is a® = ao{r, g,d, N) > such that any family A = 
{Ai, . . . , An} of covariant operators, the boundary values of the N -point Green function 
Ga defined in |2J) are real analytic in domains {E G M N ; \Ei — Ej\ > a \X\}. 

The proof will use the random walk expansion proposed by Frohlich and Spencer ^2] and 
used by Constantinescu, Frohlich, and Spencer ^U] for the DOS, which applies in the 
perturbative domain |A| <C 1 of strong disorder. As an immediate corollary we recover 
the analyticity result of Bovier, Campanino, Klein, and Perez [7]. improving the result of 
Constantinescu, Frohlich, and Spencer [10], for the case N — 1. More importantly, we 
obtain the first results on the existence of a density for the current-current correlation 
function. 

Corollary 1 For the Anderson model, under the assumption of Theorem^ for any e > 0, 
there exists a A rj£ > so that for |A| < \ r , t , the DOS is analytic in a strip of width r— e > 
about the real axis. Furthermore, there is a constant < 02 < 00 so that the current- 
current correlation functions are real analytic in (E,E') on M?\{(E,E') | \E — E'\ < 
a 2 \X\}. 

It is not expected that the current-current correlation function m(E, E') be analytic near 
the diagonal. In fact, a perturbative approach [TH HO] predicts m(£a,i?2) — C\E\ — 
E2\ 2 \n. d+l \E\ — £2 1 , and a recent result by Klein, Lenoble, and Miiller ^2] proved an 
upper bound on the AC conductivity in the spirit of Mott's formula. This has the form 
&e{v) < Cqv 2 log(-), as v — > 0, for energies E in the localization regime at strong disorder. 
The averaged conductivity OEiy) is defined as z/ _1 S s ([0, v\) for a well-defined conductivity 
measure (see J3])- If the current-current correlation function has a density miEx^ E 2 ) 
near the diagonal E\ m E 2 , then this result implies m(Ei, E 2 ) < C\E\— E 2 \ 2 \n d+2 \Ei— E 2 \ 
in the localization regime at strong disorder. The assumption that such a density exists 
near the diagonal is still unproven. We also have reason to believe that this upper bound 
might be supplemented by a similar lower bound, namely 

Problem 1 Prove that the current- current correlation function is given by a density that 
vanishes at E 1 = E 2 like \E± — E 2 \ 2 hi a \E X — E 2 \, for some constant a > depending on 
the dimension d. □ 



7 



We mention that this behavior is not consistent with the behavior predicted by the 
Mott variable range hopping conductivity argument fSl EH] , which should, in principle, 
imply the existence of an essential singularity at coincident energies. If so, it casts some 
doubt on the ability of the one-particle Anderson model to account for the properties 
of semiconductors at very low temperatures. We intend to discuss this behavior for the 
2-point, and the general iV-point correlation functions, in the strong localization regime 
in a companion article [Hj. 



3 Random Walk Expansion of the Resolvent 

The method used in these notes is the one proposed in the early days of the Anderson 
model by Frohlich and Spencer ^2] and used to treat the DOS for Gaussian distributions 
by Constantinescu, Frohlich, and Spencer JU]. It is a simple perturbation expansion in 
the small parameter A appearing in the definition (JIJ. From the definition of Ga{z) in 
flBJ), it follows that 

f 1 1 

G A (z)= / <flP(o>) (0|— \x 1 ){x 1 \A x ^\y 1 ) ■ ■ ■ (y N -i\— \x N )(x N \A NtUJ \0), (9) 

Jn — Z\ n^ — z N 

where repeated coordinates are summed over all Z d . We first show how a general Ap- 
point function Ga{z) can be approximated by an iV-point function constructed from 
simpler covariant operators A that we call r -monomials. These are effectively finite-range 
operators with analytic coefficients. 

Definition 1 Given r > 0, an element A e L°°{%>) is called an r-monomial if there is a 
family bi, ■ ■ ■ , of bounded holomorphic complex-valued functions in the strip B r = {z e 
C; \Qz\ < r}, vanishing at infinity, and a finite set of points Ui, ■ ■ ■ ,ul e Z d , such that 
the matrix elements of A^ satisfy 

L 

(0\A UJ \x) = l[b J (V(u J ))5 XtU0 . 

i=l 

A covariant operator A is an r-polynomial if it is a finite sum of r -monomials. 

Examples of covariant operators A^ include the Laplacian A = A, the velocity operators 
Vj = i[H u ,Xj], and inverses of even polynomials in the random potential A = with 
positive coefficients. These are not r-monomials since the coefficients don't vanish at 
infinity. However, they can be well-approximated by finite linear combinations of r- 
monomials as the next proposition shows. We call a covariant operator A^ a finite-range 
covariant operator if there is a finite number R > such that if x, y G Z d with \x — y\ > R, 
then (xlA^ly) = almost surely. 
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Proposition 1 Given a family A of N elements in L°°(Tf), the N -point function Ga 
can be approximated uniformly on any compact subset of (C\M) 7V by a sequence of linear 
combinations of N -point functions involving only r -monomials. 

The proof involves the following steps. 

Lemma 1 Let A be an element of the L°°{%>). Then, 

(i) its matrix elements can be written as (x\A^\y) = a(T~ x u,y — x) where, for eachu G Z d ; 
the map a u : uo G Q i— > a(cu,u) belongs to L°°(0, P). 

(ii) the sum |«J 2 converges in L°°(Q,¥). 

Proof: The covariance condition applied to A implies (x\A u} \y) = (x — s\A T - SuJ \y — s) 
P-almost surely for all s G 7L d . If we define a u (u) = a(u,u) = (0|v4Ju), then choosing 
s = x leads to (x\A u \y) = a(T~ x ui,y — x). Since uj G Q i— ► A w is measurable, so are each 
of the maps a u . Since A is bounded, it follows that ||o u ||loo < ||A||. Hence a u G L°°(Q,¥). 
In the same way, one shows that (0|A w A* |0) = J2 u ez d ku^)! 2 ^ ll^ll 2 ; showing that this 
sum actually converges in L°°(Q, P). □ 

Lemma 2 Any element A G L°°(7ip) can be weakly approximated by a sequence of r- 
polynomials. 

Proof: Thanks to Lemma ^ii), A can be uniformly approximated by a finite range 
operator, namely given e > 0, there is iV G N such that ^2^ >N \a u \ 2 < e 2 . Hence 

setting = a u , if \u\ < N, and = 0, otherwise, we obtain a finite range operator 
Ajv G L°°(7p) such that \\A — A^\\ < e. On the other hand, since Q can be taken as the 
Cartesian product M. z , any element of L°°(f2,P) can be weakly approximated (for the 
weak-* topology) by a continuous cylindrical function. A function / : O i — C is called 
cylindrical if there is L G N*, a continuous function F G Co(M L ), vanishing at infinity, and 
a finite subset {xi, ■ ■ ■ ,xi} C Z d such that f(u) = F(V u (xx), ■ ■ ■ , V u (xl)). Since C (M L ) 
is the uniform closure of the algebraic tensor product Co(M)® L such a function F can be 
uniformly approximated by a finite sum of functions of the form fiiV^Xi)) ■ ■ ■ 
where f\ G Cq(M). At last, the space 7i r , made of functions in Co(M) that can be continued 
as holomorphic functions on the strip B r vanishing at infinity, is dense in Co(M). □ 

We let H, be the space of functions on (C \ M) N that are holomorphic and bounded at 
infinity, and endowed with the topology of uniform convergence on compact sets. 

Lemma 3 The map A = (Ax,-- - ,A N ) G L°°(T P ) xN \— > Ga G Tt is multilinear and 
weak-* continuous. 

Proof: This is a consequence of the GNS-representation theorem and the property that 
on a von Neumann algebra the weak and the strong topology coincide on bounded sets. 
□ 
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Proof of Proposition [U It is a consequence of the previous Lemmas ^ El and El d 

We now return to the expansion (J2J) • Since the kinetic term in (jQ is given by the discrete 
Laplacian (without diagonal term), the Green function can then be expanded in formal 
power series in A. This gives 



oo 



3=0 



Hn 



'V 



3 

\x) . 



We next use the fact that the Laplacian H couples only nearest-neighbor terms to obtain 
a path expansion of matrix elements of the resolvent. We need some notation. We denote 
by 7 a path from y to x, namely it is a sequence (x = y, x%, ■ ■ ■ ,x n _i,x n = x) where 
Xk E Hj d for all fc's and \xk — arjfc-il = 1, for 1 < k < n. Note that all the points need not 
be distinct. The points x and x n are called the initial and the final points of 7 and will 
be denoted by <9o7 and <9i7, respectively. We call n = \j\ the length of 7. We denote by 
^(7) = { x k j < k < n} the family of (distinct) vertices of 7 4 . It will be convenient to 
denote by #7 the cardinality of V(7)so that ^7 < n + 1. The path expansion takes the 
form: 

n 

(^i^E^n^^ do) 

Given u E Z d , let n 7 (u) be the number of G [0, n] such that = u, so that n + 1 = 
6Zd n 7 (u). Hence (llOj) can be written as 

This formula will be used to expand Ga(jl) given in ©. Since Ga(z) contains iV such 
Green functions, this expansion will require N paths, namely a family r of N paths, 
T = (71, • ■ ■ ,7at). For such a family, the following notation will be used: |T| = J2iLi | Ti 1 5 
V(r) = U? =1 V(ji), #r is the cardinality of V(r), while if u E Z d , n r (u) = T,?=i n jA u )- 
The initial and final points of T are defined by d T = 9 7i an d diT = <?i7at, respectively. 
An A^-path is A compatible if (9i7i| A^l^oTi+i) 7^ 0, for i E [l,N], with the convention 
7at + i = 71 . Let then 7(A) be the set of A-compatible A^-paths with <9or = diT = 0. 
Thanks to the Appoint correlation can be written as 



AT JV 

0^)= jFiu) 2 (-A)' r ' n n (m _ - nw^iw • 



4 If the path 7 passes through the same vertex more than once, V(7) has less than n + 1 elements in 
general. 
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Let the operators Ai )U) be all r-monomials. It follows from Definition ^ and Lemma Q that 
their matrix elements factorize according to 

(x\A i>UJ \y) = S y ^ x>Ui \ [ a i>u - x (V(u)) , 

for some Ui G Z d and where a i u = 1 for all but a finite number of indices (i, u). It follows 
that an A^-path is A-compatible if and only if 

r e 7(A) ^ d oli+l - d xli = Ui Vie [1, N] . (12) 

Hence there is no need to insert the product of the Kronecker symbols associated with 
the matrix elements of the A^'s. Then, since the probability measure P factorizes with 
respect to the sites, the previous formula becomes: 



JV 

G A {z) = ^ ~ 

re7(A) uei d " ~°° i=i 



As in (J21I) of Section El if n = (ni, • • • ,njv) G N w and z = (zi, • • • , z N ) G (C \ E)^, it is 
convenient to set 



/+oo N 
^ dv h ( v ) II _ Zi y l+ i 



(13) 



Using the notation n T (u) = (n 7l (it) , • • • ,n lN (u)) and / = (1, 1, • • • ,1), the previous 
formal A^-path expansion becomes 

Ga{z) = (" A )' r| II JnAu)-f{9T,u\z), (14) 

where 

N 

9r,u( v ) = 9(v)\\ai^ dll Xv) ■ 

i=l 

As a final result in this section, we estimate the number of A-compatible A^-paths of a 
given length. 

Lemma 4 Let A be a family of N r-monomials. Then the number of A-compatible N- 
paths of total length n is bounded from above by 

#{r G 7(A); |r| =n} < (2d) n . 



11 



Proof: The initial point of T is fixed at x — 0. There are 2d ways of choosing a neighbor 
of 0. This gives the first vertex xi of 71. If xi,X2, ■ ■ ■ ,Xj have been chosen, there are 
again 2d neighbors of Xj giving 2d admissible choices for Xj + \. Hence there are exactly 
(2d)' 71 ' ways of choosing 71. Assume 71, • • • , 7fc have been chosen. Then, the final point 
<9i7fc is fixed so that the ^-compatibility fixes the initial point <9 7fc+i unambiguously due 
to (|12j). The same argument shows that the number of possible choices for 7& +1 is at most 
(2d)l 7fe+1 l. This leads to the result by recursion on k G [l,iV]. □ 



4 Estimates on Cauchy-type Integrals 

In order to prove analyticity of the correlation functions Ga{z), as expressed in (JHJ), we 
first need to analyze the Cauchy-type integrals Jn(h; z) defined in (fT3|) . We begin with 
estimates on the simplest form of these integrals for which n and z depend on one variable 
only. We will treat the general case in section 6 after we apply the results of this section 
to the DOS. The first result on the behavior of the density g is the following 



Lemma 5 Let g G l~L r and let be its n th derivative, then for any < p < r, we have 

1 I, 1 n /, n \ 7^!||n 

SUp \g(z)\ < \\g\\ r , \\g K '\\r-p < — \\g\\r- 

z:\Qz\<r-p p n 



Proof: Let z G £> r -p = {z E C \ \$$z\ < r — p}, and let 7 denote a path contained in the 
strip B r homotopic to the circle centered at z of radius p. Thanks to the Cauchy formula 



9{z) 



dz' g(z') 



2tir (z' — z) 

Taking 7 as the union of the lines 7± = {z + u ± ip ; u G R} this gives 



\9W 



du (g(u + z — ip) g{u + z + ip) 



1 ,, 
< — 9 

TVp 



2tTT I U — ip U + ip 

Using again the Cauchy formula, with now 7 being the circle centered at z of radius p, 
gives 



Integrating the absolute values of both sides over a line parallel to the real axis gives the 
result. □ 



Let I n (g] z) be defined by 
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which is convergent for ^sz ^ 0. We have the following identities. 
Lemma 6 If g G H r , then for any z G C with Ssz ^ 0, 

/»(</;*) = —his-,*) = ^W n) ;z). 

Proof: The first identity is a direct consequence of the definition. For the second 
equality, we note that 

dv 1 1 



(v — z) n+1 n (v — z) n 
Since g G H r , it follows from an integration by parts that 



n — z) n n .l-oc K v ~ z ) 



The first term vanishes while the second gives I n (g;z) = I n ^i(g^;z)/n. The formula 
follows by recursion. □ 

Lemma 7 If g eH r , then 

lim I (g,E±*) = r du 9(E + u)-g(E-u) ± 

Proof: Thanks to Cauchy's formula, this limit can be computed by using a deformed 
path 7 avoiding the point z = E. A possible choice, for positive imaginary part, is 
7 = 7_ U 70 U 7+, where 7_ = (— oo, e], 70 = {ee 10 ; —ir < 9 < 0}, and 7+ = [e, +00). This 
gives a decomposition of Io(g; E + lO) into three integrals J_ + Jo + /+. After the change 
of variable v = E ± u, the contributions of j± are given by 

f°° du /rn , f°°du ^ 

Setting v — E + ee ld gives 

/o 

Since g(i? + w) — g(-E — w) vanishes like 0{u) for tt — > 0, the sum I + + J_ converge as 
e — > giving the result. □ 

Lemma 8 For r > 0, if g G 7i r , £/ien 

sup |/ (o; £± z0)j < ( (- + 2)1 + - + 1 ) ||o|| r . 
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Proof: Using the formula given in Lemma the integral (j!6|) over u decomposes into 
Jq du(-) + du(-). For the first integral, we write g(E + u)—g(E — u) = J^ u g'(E + x) dx, 
and integrate by parts in the variable u, to obtain, 

dx g w (E + x) 



1 du g(E + u)-g(E-u) 



Another integration by parts, gives 



1 du 
(i u 



f 

Jo 



du 



g{E + u)-g{E-u) 



u 



-- [ dx In {1 / x) (g {1) {E + x) + g (1 \E -x)) . 
Jo 

( g 0)(E+l)+gP>{E-l)) 

- / dx(g {2) (E + x)-g {2) (E-x))( x — xlnx)(17) 
Jo 



To bound this, we use Lemma and obtain, 



1 



\ 9 {1) (E')\<- dlff (1) ||r-p< 



ME' E R, < p < r . 



7r(r — p) np(r — p) 

The choice p = r/2 gives the optimal bound on the right side, so that 

sup \g {1 \E')\<— 2 \\g\\ r . 
E'm tit 1 

Since < (x — xlnx) < 1, for < x < 1, the second term in (fTTj) is dominated by Hg^H? 
for all < r' < r, so that, using Lemma El again, 



■ fl (£ + M )- 8 (g- M ) 



ti 



<(- + 2) 3 



7r 



The second part is simply dominated by 

g{E + u)-g{E-u) 



du 



u 



< 



Thanks to Lemma El again, the last term on the right in (|T6l) is bounded by 

1 



Overall, this gives 



\mg{E)\ < 



\U<,:E±,())\ < ((1 + 2)1 + 1 + 1 

7r r z r 



□ 
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5 Analyticity of the Density of States 



In order to illustrate our technique in the simplest setting, we prove the analyticity of the 
density of states (DOS) in the strong disorder regime. This is a new and different proof 
of a result in [Zj obtained that was obtained using the supersymmetric replica trick for 
the Green's function. It is an improvement of the result of Constantinescu, Frohlich, and 
Spencer [10J. The case N — 1 requires only the Cauchy-type integral estimates of section 
HI The DOS exists as a function in Lj oc (M.) for lattice models since the integrated density 
of states is globally Lipschitz continuous. The main result of this section is the following 
theorem. 

Theorem 2 Let be the Anderson Hamiltonian defined in (QJ) with a distribution g G 
7i r , for some r > 0. Then, for any < e < r, there is A rjf > such that for all 
< A < A rj0 the DOS for is analytic in a strip of width r — e > about the real axis. 

In order to prove Theorem we recall that the DOS p(E) is given by 



using the notation of (jS)), and Puj(E) is the spectral function for introduced in (J2J). 
We consider functions the G\(E ± ie) = Gf(E) appearing in (|T£|) . We prove that each 
function is real analytic on R and that each has a holomorphic continuation to the lower, 
respectively, upper, half complex plane in a strip of size r — e, provided |A| is small enough. 
Using the random walk expansion described in section 2, we find 



p{E) 




E{(0|3G(£ + i0)|0)} 

lim 1 {Gi(E + ie) - G L (E - ie)} , 



(18) 




(19) 



7:0^0 



where the Cauchy-type integral I ni {u){g] A is given by 




as defined in fT5|) . With reference to Section® we note that 



J±(g;E)=\im±-I (gW;E±ie). 

e-»0 n\ 



(20) 
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Lemma 9 The function J^(g', E), defined in A2U\) . is analytic on \$sE\ < r and its deriva- 
tives have the form 

T\m Jn ^ E) ~ -^r l0 [j^TW. ' E ) ' (21) 



where 



I^(h; E) =\im I (h;E±ie). 



Proof. The formula (|2"TJl follows from Lemma El and Lemma [7| The analyticity follows 
from the hypotheses on g and the explicit formula in Lemma □ 

In order to prove real analyticity, we prove that the sum of terms on the left side of (|2"TJl 
converges uniformly. For this it is sufficient to estimate the sum 



5 

1=0 



f22) 



for < 5 < r. Since the diameter 5 will be shown to be independent of the energy, it 
follows that the sum converges to an analytic function for \5sE\ < 5. 



Lemma 10 For < 5 < r, the following estimate holds 

e 2 C (n + 3) 2 
4 (r - 5) n+3 



± ^re 2 C (n + 3) 2 

J n [g; •) * < —r-jz — s=5 \\9\ 



r j 



withC = (8/tt) + 2 + r + r 2 . 

Proof. Lemma Lemma |H1 and then Lemma El lead to 



I T ± ( 9 {l+n) \„ < 1 gjjg 
10 V> + /)!' V l|o °- (r-r 



r — rij 2 ' 

for any < r x < r. Summing over I and using formula (|35|) leads to a convergent sum 
provided 5 <T\ < r such that 

ll^0ll.<^pr^. (23) 



Maximizing the denominator on the right side of ()23|) over ri, gives the result. □ 

Proof of Theorem 2. Anticipating Section [71 Lemma ITB1 implies that the norm || ■ \\$, 
defined in eq. (J22J), is multiplicative. Thus, from the path expansion (JT9"|). it follows that 

iiGt(oii*< e w w nii<(«)-iii'- 

16 



The bound (n/2 + l) 2 < e n implies H^^n^u) /2 + 1) < e^L Definition (J22J) and 
Lemma ITU1 lead to 



igJ(oii*< E i a i 

7:0^0 



hi 



re 2 C||(7|| r l |V(7)l 



(r - 5) 2 

Using the inequality | V(-y) I < M + 1 this leads to 



,l7l 



(r-5)\ 



\\Gf(-)\\ 5 < 



re 2 C\\g\\ 

(r - sy 



E 



|A|re 3 C||g| 
(r - 5) 3 



7:0^0 

Finally, from Lemma EJ the number of paths of length n is bounded by (2d) n so the sum 
over paths can be changed into a sum over n to get 



|Gt(-)ll*< 



re 2 C||g|| r 
(r - 6f 



E 

n>0 L 



2d|A|re 3 C||#| 



(r - 5f 



So for any small e > 0, if A r>e = e 3 /(2<ie 3 rC , ||g|| r ), then for any |A| < A rj£ , the DOS is 
analytic in the strip \$sE\ < r — e. □ 



6 Estimates on TV-point Cauchy-type Integrals 

In this section, the general iV-point Cauchy-type integrals appearing in (|13|) - (jl4j) is con- 
sidered. Results similar to those of section 3 for the case of N = 1 will be obtained. Let 
n = (ni, • • • , tin) € N N , and let z = (zi, • • • , Zjsi) 6 (C \ M) N . The following integral will 
be considered 

Jn{g;z) = r°°dv- . — ^ — — , (24) 

with the following usual convention for multi-indices 

N N N Qn k 

n\ = Y[n k \, \n\ = Y, n k, df = JJ 

k=l k=l k=l k 

Lemma 11 The following formula hold 

■Ug;z) = - d^J^z), (25) 



and 

n 

Zi - z 



- 1 

Jo(9\ £) = E II —. *) (26) 



i=l 
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Proof: The first formula ()25j) is obtained by using repetitively the identity 

1 1 d n 1 



(v — z) n+1 n\ dz n v — z 
The formula (|26|) is obtained from the following identity, valid if z^ ^ z\ for k ^ I 



(27) 



(V - Z X ) ■ ■ ■ (V - Z N ) ^ 1 } Zi - Zj V- Zi 

v 1=1 J 

□ 

Lemma 12 The following formula holds 

U9-,z) = t E n ( " irj( 7 J , + nj)! 7 A^^)- (29) 

' - LJ - mAnA [Zi — z 7 -) m J+ n J +i mj! 

i=l m]\m\=ni j J J \ 

Proof: Applying (|2"K|) to (f27)j) gives 



1=1 vj:^« 

The following multiple variable generalization of the Leibnitz formula will be used 

\£-m*)-m*) = e tt^t?-^ (3°) 

r! cb r w v 7 ^ ll m k \ dz m * v 7 

m;|m|=r k=l 

Applying it to the previous formula, together with (|2"7|). leads to the result. □ 

The analyticity properties of these functions and their boundary values is now investigated 
in a manner similar to Lemma For q_ £ {+, — } N , let J^{g~, E) be defined by 

Jni9\E) = HmJn(g;E 1 + ia 1 ex,---,E N +ia N e N ). 
Then, as a consequence of (J2T?|) 

Lemma 13 If a £ {+, — j 7 ^, £/ie function E £ E i— »■ J^(g; E) £ C analytic away from 
coincident points. Moreover, away from coincident points, 



i=l m;|?n|=(i+ni 



J= 

nj + L 



-l)^( m]+ l ]+nj )\ 1 
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where 



I^ih; Ei) = ]hal (h; E { + ia^). 

Proof: Using equations (}2"o'j) and ()29|) leads directly to this formula. The real analyticity 
comes from the use of the contour of integration used in the proof of Lemma provided 
each semicircle has radius satisfying < 6j < r. □ 

Let Aat be the usual (N — l)-simplex, namely, the set of s — (si, • • • , sn) G [0, 1]^ such 

that \s\ — s\ H h Sat = 1- Let d N_1 s denote the measure ds\ ■ ■ ■ dsN-i defined on An, 

and set s- = rifcLi s k k ■ The following result generalizes Lemma [7| and eq. HdI 

Proposition 2 The function E_ G i— > E_) G C zs analytic away from coincident 
points and 

fc=l m,|mj=rifc j^fc 

where E g i— ► J = J™ 9 (g;E) E C is analytic everywhere and given by 

J= / r A , g " ,+M -'(£ t '^+«)-^ +M -"(£ t '^-), (32) 

In addition, if G\ = ■ ■ ■ = ctn — ±1, i/ie function E G R^ i— > J^(g;E) & C is analytic 
everywhere and given by 

ai = ■ ■ ■ = a N = ±i => Jf (s; £) = J^G/; ± ^K(<?; , (33) 

i^fotf) = / a ^ (§) <? (JV+Ishl) • ( 34 ) 

□ 

Remark 1 The previous proposition shows that admits a polar singularity at Ef. = Ej 
if and only if ^ <jj . 

The first step in the proof is a simple lemma relating integration over an N — 1-simplex 
to the product of the singular terms in the integrand in ([24)1 . Given a set of points 
{z\, ■ ■ ■ , zn}, their convex hull will be denoted by Conv{zi, • • ■ , zn}- 
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Lemma 14 Let n G and z G (C \ R)^. If v ^ Convjzi, • • • , zjv}, t/ie following 
formula holds 

n 1 = (^+n- 1 )' / ^-1, g 

11 („ - 2fc )n fc+ l a | 7 Ajv - (u _ s ^)iV +N ' 

Proof: Both sides of this formula are defined and holomorphic in v in the domain 
{v G C ; v ^ Convjzx, • • ■ , Zjv}}. Thus, using the unique analytic continuation theorem, 
it is sufficient to prove it for 3ft(t> — z fe ) > 0. Setting a k = v — z k , the identity 

i poo ±n k 

— — = / dt k ^- e~ tka \ 

Oj k Jo 'Ik- 

valid for $ta k > 0, implies 

N 

11 ( afc )n*+i y R? - n! 
The following change of variables (t±, • ■ ■ , tjy) i— > (A, si, • • • , sjy-i) will be useful 

AT-1 

= As fc , (1<A;<JV) sjv = 1 - ^s fc , A > 0. 

fc=i 

In particular, we have s k > V7c, and Sj + • • • + Siv-i < 1. Furthermore, the volume 
elements transform as 

dtiA---Adt N = (-l) N ~ 1 \ N ~ 1 d\Ads 1 A---Ads N - 1 , 

so that 



^ 1 /■ n /-oo 

Hr^TT = / ds 1 ---ds N . 1 ^ / A^^dAe^^ 



(N+\n\ - 1)! 

asi ■ • • dsN-i — 7 



Replacing by t> — Zfc gives the result. □ 
Proof of Proposition 

1. The analyticity claim in the first part of the proposition follows once we have proved 
the representations (J31|) and (|32|). If Zi, ■ ■ ■ , Zjj are all on the same side of the real axis, 
then the convex hull is also contained in the same half plane. Thanks to Lemma El it 
follows that 
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From Lemma [7[ equations (jH2J), fjHHj) . and (}3"Tj) follow immediately. 

2. To derive (}3*T]) from in the case ai = ■ ■ ■ = a N = ±1, we apply Lemma ITU for 
ni = • • • = njv = 1, together with (J28|) . to obtain 



Ja n -(z-I2 k SkE k ) N {z-Ek)...{z-E N ) j^z-E k ^j k E k -E - 

Multiplying both sides by g(z) and integrating over a Jordan path surrounding each E k 
once and contained in the holomorphy domain of g, leads to 

Applying the differential operator (l/n\) d\ on both sides gives 



An application of Leibnitz rule (see (J30J)) gives equation (|3~Tj) for o\ = ■ ■ ■ = = 1. 

3. If a = (oi, • • • , ctjv) ^ {+1) — 1}^ is arbitrary, let A(q) C [1, N] be the set of indices k 
such that crfc = — 1. Whenever 2 = • ■ • , zn) £ is such that a^z^ > for all fc's, 
the integration contour defining Jn(g;z) cannot be deformed to have all z k s on the same 
side. Actually the 2%'s such that k G A(cr) are below the integration path while the other 
2%'s are above. Deforming the contour to get all the z k s above can be done to the price 
of adding the corresponding residues, leading to 



JS-(E) = J+(£)-2z7r V Res z=Eh - — — ^ '-. . - 

nK—J nK—J * k (z - E^^ 1 ■ ■ ■ (z - E N ) n " +1 

= J~(E) + in > (ai. — l)Res z -R, — — ; -— — . 

= ( z - #i) ni+1 ■ ■ ■ - ^) niV+1 

The residues are given by 
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y.+k 



Use of the Leibnitz rule again together with the previous formula for Jn(g',E) gives (J3T|) 
in the general case. □ 



7 Real Analyticity 

We now use the estimates on the Cauchy integrals obtained in Section |S] in order to 
prove that the iV-point Cauchy-type integrals, for N > 2, are real analytic as functions 
of £ 6 N , on the domain for which Ej ^ Ei, for i j. To this end, we prove that 
the functions J^(g;E), defined in section 6, have uniformly convergent power series on 
appropriate domains. It is convenient to do this using certain Banach spaces of real 
analytic functions. Let 5 > and let U be an open subset of M. N . For a continuous 
complex valued function / on U let ||/||w be the sup- norm ||/||« = sup £gW \f(E)\. Then 
let Ss(U) be the space of smooth functions F : E_eU i— > C vanishing at infinity, together 
with all derivatives, such that 

S\l\ 

\\F\\s,u = Yl 7T H^U" < 00 • 

l€N N ~ 

Endowed with the norm || ■ \\gu the space EgiU) is a Banach space. Moreover, it is a 
Banach *-algebra if endowed with the pointwise multiplication and complex conjugacy 
since 

Lemma 15 Let F and G be two elements of £§(JA). Then both F* and F ■ G belong to 
£s(U) and 

\\F ■ G\\$ t u < IkIUm ll^nkw ) ll-F*lkw — ll-nkw 

Proof: The only non trivial property is the first inequality. Using the Leibnitz rule (see 
(UP) leads to 

id l -(F-G) = V i- d^F i- d-G . 
/! ^ ml nl 

m,n;m+n=l 

From this the inequality follows immediately since 5 1— + -' = <5'— ' 5'-'. □ 
When we take IA = M, N , we write || • ||^ in place of || • \\s,rn- 
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Proposition 3 Let g be an element ofH r , for some r > 0, let C = 8/ir + 2 + r + r 2 , 
and consider N > 2. For any choice of S > and A > such that < 5 < A/2 and 
0<A — 5 < r/2, the boundary values of J%(g) satisfy the following estimate 

\\ J n\.9)\\S,V A < 



r \A-5 / 
in the domain V A = {E G R N ; \E { - Ej\ > A}. 

The proof will proceed in several steps. First the following identity will be used often 

(1 = f(i±^ ul . (35) 

1=0 

The following combinatorial estimate is necessary. 
Lemma 16 If n e N L and if r GN, i/ien 

E (m + «,)! (r + \n\ + L - 1)! 
m! n! ~~ r! (Inl +L - 1)! 

Proof: If C^(n) denotes the left side of (jHEJ), its definition and the identity (|33|) give 

E = n E - (13^- 

r=0 fc = lm fe =0 V 7 

Expanding the right side in formal power series in X by using (|33|) again gives directly 
the result. □ 

Proof of Proposition OH Using Lemma IT31 and assuming \Ei — Ej\ > A gives 

A ( (rrij + lj + Uj)\ ( 5\ lj 1 



»,EE E n^^({) 



»=1 ZeN-^ m;|m|=2i+rii 

( / 7._L n . v )l , / ("»i) 



A 7 A m J +n J +1 



oo • 



Using the distributive property of the product with respect to addition for the variable lj 
gives 

X x x> (nij + lj + nj)\ f 5 V 



m)\u,v> £ EE E n E 



r-l I,=() nr. m_ ^/,+», I./:--' TTlj . lj . Tlj . 



Jtj + ^ Sk \\l ( 9 — 
IA nA. \ mA 
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We use Lemmas 03 and |H1 to estimate the L°°-norm in (|37|) for any < r j < r by 

\ ffij! y (r — rjj^ TOj! n 

Using to sum over /j, together with estimate (|3*SJ). leads to 

kmii^ < EE E (n (m ' +K))! 



m,! n-,! (A - <5) m J +™ J + 1 

j=l i i= m-,\rn\=h+ni 

x- 



2 llJ/llr > 



for < rj < r. We define p by |m| = m« +p and |n|j = |n| — fii. With this definition and 
Lemma EE we perform the restricted sum over m G N and obtain 

\\r(a)\\ < f 1 V + r< V 

IKn^||5,o A _ Z^( A _ 5 ) Ni+ iv-i^ LL p!(|n|i + JV-2)! 



X 



WA- <y)Hi+' 

1 v 7 Z i= \ p=0 



raj r ™* +ii (r — Tj) 2 
Using (jHHj) again to sum over p, if < A — 5, leads to 



\\JnK9)Uv A _ ^ {r _ ri)2r n i{A _ 5 _ ri y nU+ N-l2^ y n .j 

Finally, using again (J35|) to sum over Zj, if 5 < < r, leads to 



\\ T gj a \\\ < V Cr h\\r 

|| n\y)\\5,-D A _ ( r _ r .)2( r ._ 5 )n i+ l( A _ 5 _ r .)]„| i+ iV-l " 



These estimates are satisfied on the domain provided < 5 < < r, and (A — 5) > 0. 
In particular, if we take A > 25. Recalling that we treat here N > 2, if A — 5 < r/2, it 
follows that Ti < r/2, so that r/(r — r^) 2 < 4/r. Then, minimizing over in each term 
on the right side gives 



leading to 



K + 1)(A- 


5) + (\n\ 


i + N-l)S 




n\ 


+ N 
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n ■ 4- 1 

r,-5=(A-5)^-, A-5-n = {A-8) 

\n\ + JS 

If l,m are positive integers, ((/ + m)/l) [ = (1 + m/l) < e m , so that 

l + m\ 1 ( I + vii 



\n\ 


i + iV-1 




n 





I I \ m 



< e l+m , 



we obtain 



4A/T7IMI / p xlsl+JV 



r V A -<V 2 
This proves Proposition 3. □ 



8 Proof of Theorem U 

It is now possible to finish the proof of Theorem ^ proving the main theorem on the real 
analyticity of the iV-point correlation functions Ga{E_), away from a small neighborhood 
of the coincident points E^ = Ej. We recall the A-compatible path expansion in (|14|). Our 
strategy is to write Ga(E_) as a sum over n G N, the length of the compatible paths. We 
then prove that Ga{E_) is m the Banach space Ss(T>^) using the estimates in Proposition 3. 
We first prove the theorem for the case when the covariant operators Ai are r-monomials 
which implies the result when the covariant operators Ai are r-polynomials. Starting from 
(114}) . it is sufficient to prove that the formal path expansion converges in Ss(T>^), with the 
domain C R , for convenient values of 5 and A. Using Proposition El the Appoint 
correlation G\ is estimated as follows 



\\G%\ & ^ < l A l' r| II ^illflr,- 

re9(A) uev(r) 

Clearly, we have that \n r (u) ~ f\ + N = \llr( u )\- This is bounded above by the length of 
the path |T|. On the other hand, if A = max itU ||a» w ||, then A appears at most AT-times 
in the right hand side in the definition of gr jU , reflecting the fact that G a is homogeneous 
of degree N in A. As a consequence, we get the bound 

II NC^uWr < A N H NdWgWr. 
nev(r) «ev(r) 

On the other hand, since g defines a probability on R, \\g\\ r > 1. Moreover, C\ = 
4(8/tt + 2)/r + 4 + 4r > 4, so that iVCi||p|| r > 1. Hence, using #r < |r|, and the 
Lemma EJ this inequality becomes 



A-5 



\RrV u )-J |-t- JV 



r 
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|r| /„\\\Arri lull \ l r l 

r \ 



IGlll^ < A" £ (^L) II ^,11,11, < A" e ( 



A-5 



^ ^g-WeiAHMi.y < _ £ <oc 

n=0 ^ 



A -5 y - l-2dNC 1 e\\\\\g\\ r (A-5)- 1 

(39) 



provided 



IAI < 



A-5 



2dNC 1 e\\g\\ r 

Choosing 5 < A/2, we define a = 4diV(7ie||gr|| r . Hence, the function G £s(T>a), 
provided ao|A| < A, proving Theorem 1 for the case when A is a covariant family of 
r-monomials. We now consider a general covariant family A G L 0O (7p) x7V and the as- 
sociated correlation function Ga{z), with Qf^j- > 0. By Proposition 1, Ga(z) can be 
approximated uniformly on any compact subset of (C\M) Ar by a sequence of correlation 
functions Ga (z), where the family A n are r-polynomials. Furthermore, the coefficients 
can be chosen so that they are uniformly bounded in L°°(Q, P). This family of correlation 
functions has an analytic continuation to £>a- Estimate (J39)) show that this family of ap- 
proximating correlations functions is uniformly bounded since the bound depends on the 
coefficients. Since the family converges uniformly on any compact subset of the original 
domain (C\R) , the family converges uniformly on any compact subset of the extended 
domain in C N . The limit functions provide a continuation of Ga(z) into this extended 
domain by the identity principle for analytic functions. This proves Theorem 1. □ 
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